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1 Introduction 

A lot of problems in QCD at zero temperature, finite temperature, and in the case 
of non-equilibrium situations require the use of non-perturbative methods. Important 
examples are the mechanism of confinement, probably related to the vacuum structure 
, the investigation of hadron spectra and hadronic matter, the understanding of the 
nature of the deconfinement transition and the chiral symmetry restoration. Further- 
more, the mechanism of hadronization in relativistic heavy-ion collisions, - following 
the possible preformation of a quark-gluon plasma - or possible non-perturbative effects 
in a quark-gluon plasma above the critical temperature cannot be treated perturba- 
tively as well 0|. Besides numerically very involved lattice calculations [^ there are 
only preliminary attempts to incorporate non-perturbative effects by using the Dyson- 
Schwinger equation or variational calculations @] - .In addition, considering the 
formation, evolution, and decay of a quark-gluon plasma in relativistic heavy-ion colli- 
sions, a transport theory based on gauge covariant Wigner functions has been proposed 
by Elze et al. ||, H). It is, in fact, a many-body theory completely equivalent to the 
Heisenberg equations for the quarks and gluon fields, however, as a reformulation of 
the SU(N) gauge theory in phase space not manageable in practice. The standard ap- 
plicable limit is the semi-classical, abelian approximation which corresponds to lowest 
order perturbation theory in the high temperature approximation. A more promising 



approach may be provided by the use of equal-time Wigner functions fTO . 

In this paper we propose a different approach along the line of relativistic two-body 
correlation dynamics [Tl|-[0, which has proven to provide the genuine basis for the 



formulation of non-perturbative transport theories for baryons and mesons |T^. It is 
already known that this approach obeys all conservation laws concerning space-time 
degrees of freedom, i.e. fermion number, linear and angular momentum as well as the 



total energy |TT|. The novel phenomena in the application to quarks and gluons is the 
non-abelian dynamics of the gluon fields as well as a SU(N) internal gauge symmetry. 
The paper is organized as follows: In Section 2 we first give a brief reminder of 
the SU(N) gauge theory in the temporal gauge and introduce the necessary notations. 
Section 3 is devoted to the derivation of equations of motion and suitable truncation 
schemes for linked n-point Green functions. In Section 4 we show the conservation 
of the Gauss law in time (as a quantum expectation value) within the approximation 
scheme adopted and investigate the higher order Ward identities (as generated by com- 
mutators of Gauss operators), whereas Section 5 provides a summary and discussion 
of open problems. The explicit form of the final constraint dynamical equations is in 
part shifted to the Appendices in view of their length. 



2 SU(N) gauge theories in the temporal gauge 



Here we repeat the basic concepts and equations of SU(N) gauge theories ||T6[, which 
build the starting point for the correlation dynamics in QCD. 
We begin with the SU(N) gauge Lagrangian 

£ = -^7,^,M/ - iF;,F;,, (2.1) 

where the gauge- covariant derivative D^ is defined as 

D, = d,- igT^Al . (2.2) 

The fields A" are the gauge potentials and T" are the N x N generators of SU(N) 
following 

TriT^T^} = -6"^ (2.4) 

[T'',T^]=ir^^T', (2.5) 

whereas the field strength tensor is given by 

F;, = d,At - d.Al + gr''A\Al (2.6) 

or in terms of color electric and magnetic fields as 

E] = iFf, = -A] + id.Al + igr'^'^A^Al, (2.7) 

Bl = y^'Ft^ = ^e'^\d,A^ - d,A\ + gr''A\A^^. (2.8) 

The Euler-Lagrange equations of motion for massless quarks then read 

7^D^^ = 7,(5^ - igT'-AD^ = (2.9) 

and 

5.C + air'^K^'^u + K) = (2.10) 

for the gluon fields, where /^ is the quark-color current 

/« = i^^^T''^. (2.11) 

In order to perform actual calculations we have to fix the gauge. In view of the 
standard Hamiltonian formulation of correlation dynamics we choose a special axial 
gauge, i.e. the temporal (or Weyl) gauge [O 



Al = 0, (2.12) 

such that the conjugate field of Af is (except for a sign) the color electric field 

Iit = M = -Et- (2.13) 



In the temporal gauge we obtain the canonical quantization conditions 

[A^ir,t),I[%r',t)] = t6,, S^' 6'ir - r'), (2.14) 

{^ir,t),^\r',t)} = 6%r~r'), (2.15) 

while all other commutators or anticommutators vanish. The Hamiltonian and its 
density then is given by {x = (r, t) ) 

H = f £rn{x), (2.16) 

with 

nix) = ^{U^W} + B^B^) - gl^A'^ + ^7,^,^ (2.17) 

and the Heisenberg equations of motion for 11" read 

n," = -[n^ H] = Df^F^^ + gll (2.18) 

t 

with 

Dl' = 5'''du + gr''Al (2.19) 

while the equations of motion for A1 reduce to (cf. (2.13)) 

At = -[AlH] = Ii1. (2.20) 

Since the Euler-Lagrange equations with z/ = cannot be derived from the Heisenberg 
equations of motion, we have to deal with additional constraints on the gauge fields, 



i.e. the Gauss law 16, 17 



g^[x) = r{x) + <^\x)T''^{x) = (2.21) 

with 

r(x) = -D'^'m. (2.22) 
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In quantum physics (2.21) cannot be fulfilled as an operator equation but has to be 
imposed as constraints on the physical states \p), 

g\x) \p) = 0. (2.23) 

These equations are similar to Slavnov- Taylor-Ward identities within the temporal 
gauge, which impose relations between different Green functions |jl8| . 



The basic idea of correlation djTiamics now is the use of dynamical equations of 
motion for equal time Green functions, which is sufficient to describe the time evolution 
of the system in a local field theory. For this purpose, we consider the operators 
(j^"'"^(a;i, • ■ ■ , x„), which are products of the gluon- field operator A^ and its canonically 
conjugate momentum II'^ as well as quark field operators \E'^ and \l/, 

G'i'^-">(zi.---,x„) =G^">(/l^,n^',v|/t,vi/), (2.24) 



where {a, a} denotes a set of color and spinor indices, respectively, n is the number 
of field operators and all time arguments are taken to be equal. The equal time Green 
functions are given by the quantum expectation value within the physical state \p) 

G^°}(xi,---,x„) = (p|G'i"'">b) (2.25) 

and the equations of motion for these Green functions follow from the Heisenberg 
equations 

^GlrHxi,---,^n) = -{p\[Glr\H]\p)- (2.26) 

at I 

Within the canonical quantization conditions (2.14) and (2.15) equation (2.26) leads 
to a coupled set of equations for the various n-point functions, which for practical 
purposes has to be truncated. 

3 Equations of motion 

The derivation of dynamical equations of motion is straightforward, however, somewhat 
lengthy. We start from the Heisenberg equations of motion (2.26) for the quantum 
expectation value of the Green functions (2.25) and first only show the lowest order 
equations within the notation adopted. We use '^ = u'l with a denoting the color index 
and (3 the spinor index, respectively. For the Dirac matrices 7^, Oj = 70 7^ we use the 
greek indices a, (3 to specify the spinor components, i.e. 7^^ = (7^)0/3, cx\p = (tt*)o/3, 
and define t^"'^ by 

where h = 1, ..., (A^^ — 1) runs over all generators, while a and c run from 1 to A^. This 
gives for the 2-point quark Green function 

^|«(^)<''(^')) = 

+c^^,,[d4ul{x)4{x')) - gt'^'^{ul{x)ul{x')A\{x'))]. (3.2) 

For 1- and 2-point gluon Green functions we get 

|(^-(^)) = (n^(^)), (3.3) 

+2gr'^dyMU^)AM)k^=y) + g'r"'r'%Al{x)Ai{x)Al{x)) 
+gr''^[dyM\{y)At{x)) - dy^{A\{x)Al{y))\^=y) 
+9t''''l'.p{<{x)ul{x)), (3.4) 

j^{At{x)AUx')) = {Ut{x)Af,{x')) + (A^(x)n«;(x')), (3.5) 



j^{Ti1{x)Ai{x')) = {Yi1{x)Yii{x'))+dl{A1{x)A!i:{x')) 

-d^^d^Ml{^)4{^')) + 2gr''^dyM\{^)^Ay)Ai{x'))\^y=.) 
+9'r'j''^Al{x)At{x)Anx)At:{x')) 
+gr'^[dy,{Aliy)AK^)At'{x')) - 9,^(4(x)A^(t/)A^;(x'))](,=.) 

+9t-''lU<i^Hi^)^t'i^')) (3-6) 



and 



j^{mx)uf, ix')) = di{Atix)ut:ix'))+di,{mx)At:ix')) 
+2gr''%Ml{x)Ai{ym:{x'))\^y=^) 

+2gr''^dy,{mx)Al{x')AUy))\iy=.') 
+g'r'j''^Al{x)At{x)Anx)Ut:{x')) 
+g'r'''r''{mx)Al{x')At{x')AUx')) 

+9r'VyAAliy)AKmt'i^')) - dyAAl{x)Al{y)Ut:{x'))]^,=.) 
+gr''^[dy,{Ut{x)Al{y)AUx')) - 9,,(n^(x)4(x')A^(t/))](,=.o 
+9t''"'llfs{uUx)ulix)Ut:ix')) + gt^''^Y:^{Ut{x)u'^{x')ul{x')). (3.7) 

For the lowest order vertices one obtains 

^|«(^)<'(^')4'(?/)) = ^{K{x)ui,{x')U^{y)) 

+<[5..(«^(^)<'(^')4(?/)) - gt''^%u^p{x)ui,{x')At{x)A^{y))] 
+«''/.[54«(^)<(^')4(l/)) - i7t^'"'^«(x)«^(x')4'(^')A'(l/))], (3.8) 

^f^{K{x)ui,{x'n{y)) = 

+<[5..(^K^)<''(^')n^(y)) - gf''^u',ix)ui,{x')AUxmy))] 

+t[dl{uUx)ui,{x')A^{y)) - dy,dy,{ul{x)ui,{x')Al{y)) 

+2gr'nA<i^)<'i^')4iy)Atiz))k.=y) 

+9'r'j''^<{x)ui,{x')AUy)At{y)A'riy))] 

+^^f ^'^K(<(x)<:(xOA^>)A^(y)) - 9,^(<(x)<:(x')4'(y)A(^))](.=.) 

+^9t''Hp.{Ki^)ui'ix')4iyK'iy)) ■ (3-9) 

In general the equations of motion for n-point Green functions couple to (n+1)- 
and (n+2)-point Green functions, as can be seen e.g. from commuting (H")" with the 
cubic and quartic terms of A" in H, such that the equations of motion form an infinite 
set that has to be truncated for practical purposes. 
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A truncation scheme which has been quite successful in the nuclear physics context 
is based on the cluster expansion of Green functions in terms of connected (correlated) 
Green functions since the connected Green functions become of minor importance with 
increasing order [11]. This has been shown explicitly in the nuclear physics context in 



19| . A nonperturbative truncation scheme that leads to a resummation of loop and 
ladder diagrams in infinite order and observes the space-time conservation laws as well 
as the weak Gauss law (see next Section) is a truncation up to 4-point Green functions. 
For the case of ground state Green functions the explicit expressions for the clus- 
ter expansions can be derived from the representation of the theory in terms of the 
generating functional of full and connected Green functions, Z[f], rj, J] and W[fi, r], J] 
respectively, which are given by pO| 



Z[r], ^,J] = J dii (^, ^, A^) exp (^i J d'x [C[^, ^, A^] + ^r] + r]^ + A^ J^)) (3.10) 
and 

Z[fi,ri,J] = expiW[f],v.J]), (3.11) 

where ghosts can be omitted according to our choice of gauge. As an example we derive 
the lowest order expansions in order to present the general concept. We start with the 
cluster expansions for the time-ordered Green functions with different time arguments: 

{TAnx)At:{x')) = j;^^JJ^fMW[n,V,J])U=r,=J=o 

= (T Anx)Af,{x')), + {Anx))Mt'ix'))c ; (3.12) 

{T ulix)ui,ix')Al{y)) = I J ^ ^^^expiW[f],v,J])\,=,=j=o 

-W[r],'n,J] 



\\5r^-{x)5f^i,{x')i5J'M 
^'^^-^FTTT^t^'^''^]) 1 A^a^^^.A ^[^^^^ A 1 exp{W[n,r], J])U=r,=j=o 



(r<(x)<:(x')A^(y)).+ (T<(x)<:(x'))e(A'(2/))c, (3.13) 



where all Green functions containing an unequal number of u and u are assumed 
to vanish. The generalization to Green functions containing conjugate gluon field 
momenta is straightforward |^. The expressions for equal time ground state Green 



functions are obtained by taking the well-defined equal time limit which yields the 
appropriate operator ordering in the cluster expansions ( |3.12| , |3.13| ). For the general 



non-equilibrium case we then simply define the connected Green functions by using the 
same form of the cluster expansion as in the ground state case. 
For the lowest order Green functions one arrives at 

{Atix)At:ix')) = {Atix))Mi:i^'))c + (A^(x)A^;(x'))e (s.u) 

{U1{X)A^:{X')) = (n^(x))e(A^;(x'))c+ (n^(x)A^;(x'))e (3.15) 

(n^(x)n^;(x')) = (n^(x)),(n^;(x'))e + (n^(x)n^;(x'))e (3.i6) 

«(x)<:(x')A'(z/)) = «(x)<:(x'))e(A^(2/))c+ «(x)<:(x')A'(?/))c (3.17) 

«(x)<:(x')n,^(l/)) = {<{x)ui,{x'MU'^iy)),+ {ul{x)ui,{x'n{y)),, (3.18) 

where now the time arguments are again taken to be equal. The explicit expressions 
for the relevant cluster expansions up to sixth order are presented in Appendix A. 

In reducing the dynamics to the 4-point level we assume that all correlated terms 
of order n > 5 can be neglected, i.e. 

Q{a,a} ^ Q f^^ ^ > 5^ ^3_^g^ 

Thus the set of equations of motion becomes closed on a finite level. 

The equations of motion for the connected Green functions are obtained by inserting 
the relevant cluster expansions into the equations for the equal time Green functions 
((p.2| ) -( p^ a.s.f.). As an example we present here the resulting equations of motion 
for the connected 1-point gluon Green functions, 

f^{At{x))^={Ut{x)U (3.20) 



~dt 



n^(x)), = a,^(A(x)), - a.,a.,(A(x))e 

+2gr"^dy,[{A\{x)A'i{y)), + {A\{x))M-{y))ch=y) 

+g'r"'r''[{Al{x)Ai{x)Al{x)), 

+ (4(^))c(^^(x)A(^))c + {Al{x)Al{x))Mi{x)), 

+{A\{x)Ai{x))Ml{x))c + {A\{x))Mi{x))Ml{x))c] 

+gr'^dy,Ul{y)A'l{x)), + {Al{y))M-{x))cV.=y) 

-gr"'dyMl{^)A^k{y))c+{Al{x))Ml{y))c](.=y) 

+gt^'^^l,p{ui{x)ul{x)),, (3.21) 



and the equations of motion for the connected 2-point gluon Green functions 
d 



^^{Anx)At, (x'))e = {mx)At, (x'))e + {AtixM (x'))e, (3.22) 



|(n^(a:)A»;(xO)c = (n^(x)n^;(x'))e 

+di{At{x)Ai:ix')), - a.,a.,(A«(x)A^;(x'))e 

+2gr'^%^[{Al{x)At:{x'))MKy))c 

+ {AlixMA-iy)Ai:ix')), + (4(^)^-(?/)^-''(^'))c](.=.) 

+g'r''r'''[{Al{x)At{x)At{x)A^:{x')), 

+(1 + Z, + Z,){{Al{x))Mtix))c{A^{x)Af,{x')), 

+ {Al{x)At{x)A^:{x'))MKx))c + {Al{x)At{x))Mtix)At:ix'))c}] 

+gr%,[{Al{y)A^:{x'))MKx))c 
-gr'''dyMAlix)A^:ix'))Mliy))c 

+{Al{x))Ml{y)A'i:{x')). + (4(^)^fc(y)^-''(^'))e](.=.) 

+ge''lUu'^{xXp{x)At:{x')),, (3.23) 

|(n^(x)n«;(x'))c = (i + P..')[5^(^-(^)n^'(^'))e-5.,a..(A^(x)n^;(x'))J 

+2g{i + Vu').r''dy,[{Al{x)u:}:{x'))MKy))c 
+(4(x)),(4^(2/)n^;(x'))e + (4(^)^?(?/)n^''(^'))c](.=.) 

+/(i + p.,)r'7'='1(4(^)^'(^)4^(^)n?>'))c 
+(1 + z, + re,){(4(x)),(4(x)),(4^(x)n^;(x'))e 

+ (4(^)^f (^)n?'(^'))e(4^(^))e + (4(^)^'(^))e(4^(^)n^;(^'))e}] 

+g{l + Vu,)r'%,[{Al{y)Ut:{x'))MKx))c 

+ {Al{yMAKx)Il^:{x')), + {Al{y)At{x)U^:{x')U^,=y^ 

-g{i + v.,)r'^dyA{Al{x)ut:{x'))Mliy))c 



+{Al{x))Mliym'i^'))c + (4(^)^^(2/)n?'(x'))J(.=,) 



+g{l + Vu'r'''Y^p{u'^{xX^{x)U^, (x'))c. (3.24) 

In (|3.23|) and ( p.24|) we have introduced two kinds of permutation operators Vw and 
%d in order to achieve an unambiguous compactification. To demonstrate the action 
of these operators we look as an example at parts of the (7^-terms in ( |3.24| ) 



= g\l + Vu')r'j'''''[{AiixMAti^MAtix)nf,ix')), 

+(^t(^))c(^^(x)).(Af(x)n^;(x'))e + (An^))c(Af(x)).(4(^)nr'(^'))c] 

= 9'r''r'^[{Al{xMAtixMAt{x)Ut:{x'))^ 

+9'r''j^'''[{Ai{x'))Mt{x'Mmx)AUx'))c 

+(4(^0)c(A(^'))c(nr(^)Af(x'))c + (A,(x'))c(A^(x'))c(n^(xX(x'))e](3.25) 

Obviously Vn' and 7^^ act in the same way by interchanging the fields or conjugate 
momenta labeled by the spatial components i, i' (in case of Vn') or by the color indices 
e,d (in case of %d)- Due to their length we have shifted the remaining equations of 
motion for the connected 3- and 4-point gluon Green functions to Appendix B. 
For the connected 2- and 4-point quark Green functions we obtain 

«^.[5..(^K^)<'(^'))c - 9t'^^{{u',ix)ui,ix')Alix)), + {u',ix)ui,ix'))M-i^))c}] 

+<'/3[5x:«(^)<(a:'))c-<7t^"'^{«(x)n^(x')A(x'))e+«(x)u^(x'))c(A(^'))c}] 

(3.26) 

and 

tj^{ui,{x')ui:{x'')u';;:^^^^^^^ 

-^t^»'^{(^^(xo<:(x'o4:(/)<(yO)c(A(xO)c+(tiMxo4(yO)c(<:(x'')<"(/)4(^'))c 
-(t.^(x>j';,(i/"))c«:(:r")<,(y')A^(x'))c}] 

-(<,(xo<"(y'OA'(2/0)c(<:(x'o«K2/0)c+(<'(:^o^K2/'))c(<:(x'')<:(2/'04'(2/'))c}] 

+(<,(xo4(yo^'(^'0)c(sKx'o<:(2/'0)c-(<.(xO<"(2/'')A'(^''))c(^^(a;'')<'(i/'))c}] 
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+(<:(xo4(yOA'(/))c(<:(:^'0«^(/))c-(<:(x')«^(/))c(<"(x'')4(i/')A'(/))c}], 

(3.27) 

which represents all equations on the pure quark sector. 

The equations of motion for the connected 3-point quark-gluon vertex Green func- 
tions are given by 

^|k(^)<'(^')A(?/))c = ^«(x)<:(x')n^(t/))e 

+aJj9.,(«^(x)<:(a:')A'(t/))c - ^t^''"^{(«^(:r)<:(x')4"(x)A^(t/)), 

+ «(x)«^(x')).(4"(x')A^(l/))c+(ti^(:r)ti^(x')A'(y))c(At°(x'))c}] (3.28) 
and 

^J^{ul{x)ni{x'n{y))c = <K(ti^(x)<:(x')n,^(y))e 

-^t^«-{(«^(x)<:(x')4"(x)n,^(2/)), + (t.^(x)<:(x'))e(4°(^)n.^(y))c 

+«''/.[5<(<(^)<(^')n^(l/))c - ^t^""'={(<(a:)«^(x')At°(x')n^(?/))c 

+(ti^(x)«^(x'))c(4"(^')n.^(y))c + «(x)«^(x')n,^(2/)).(4"(x'))c}] 

+z2gr^'^d,^[{ul{x)ui,{x')A\^{y)At{z)), 

+ (^^^(^)<'(^')4"(?/))c(A^(^))c + «(x)<:(x')A^(^))c(4"(y))c](.=.) 

+ig' f'^'" f^i + r,,, + r,„e)[«(x)<:(x')4"(i/))c(4(?/)A^(i/))c 

+ (ti^(x)<:(a:0A^"(y))c(Af(y))c(A^(y))c+(ti^(x)<:(x')A^(y)A^(?/)).(At"(2/))d 
+.<7/'«^9,J«(x)<:(a:')4"(^)A^(l/))c 

+ «(x)<:(x')A^"(^))e(A(t/))c + «(x)<:(x')A(l/))c(At"(^))e](,=.) 

-z5r"^a,j«(x)<:(x')4"(y)A(^))e 

+«(x)<:(x')At°(y)),(A^(z)), + «(x)<:(x')4(^))c(4"(t/))c](,=.) 

11 



Mgt''^^ll,\{ul{x)u'^{y)ul,{y)ui{x')), - {ul{x)u^y)),{u'^{y)ui{x'))l (3.29) 

while the exphcit equations of motion for the connected 4-point mixed quark and gluon 
Green functions are shifted to Appendix B again. 

The closed set of equations ([OOl ) - (3.24) and (3.26) - (3.29) including ( pJ|) to 
(p3.12| ) will be denoted as Constraint QCD (CQCD) equations furtheron. 

4 Compatibility with Gauss law and Ward 
identities 

The aim of the present section is to show the compatibility of the CQCD equations 
with the Gauss law and the Ward identities, respectively (cf. [0). For this purpose, 
we first note that the operators g°'{x) at equal times constitute a set of local generators 
for the SU(N) gauge group because of ||16[ 

[g\v;t),g\v\t)]=ir''g'^{v,t) 6\r-r'). (4.1) 



Furthermore, the Gauss law operator commutes with the Hamiltonian |T6[ 

[g''ir,t),H] = (4.2) 

which expresses the conservation of g'^ in time 

^g^{r,t) = -[g^{r,t),H] = 0. (4.3) 

at I 

Equations (|4.2| ) and ( [4.3|) imply that within the temporal gauge the system evolved 
in time by H has a residual gauge symmetry and the Gauss operators (/"(r, t) are the 
local generators of the residual gauge symmetry. 

Now we assume the actual realization of the Gauss law as a quantum expectation 
value 

(p|^"(r,t)b) = (4.4) 

instead of (2.23). This may be considered as a weak form of the Gauss law, since it 
is less restrictive for the physical states than (2.23). For example, the perturbative 
vacuum |0) fulfills (|4^ ) but not (2.23) p3[|. Also (|4^) allows for the local propagation 
of colored objects, whereas (2.23) restricts to color-singlet objects. These consequences 



of ( [4.4| ) are exactly what we aim at for describing properties of the quark-gluon plasma. 
Furthermore, since all g°'{x) are conserved quantities and the field equations of mo- 
tion are generated by H, the weak Gauss law is conserved identically in time provided 
that it is fulfilled initially. This is just what we have been aiming at: The Hamiltonian 
dynamics are compatible with the conservation of the weak Gauss law, and the prob- 
lem is shifted to an initial value problem. This statement still holds within 4-point 
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correlation dynamics, since the Gauss law operator contains two-point operators at 
most. 

Whereas the weak Gauss law is a constraint imposed on the lowest order Green 
functions of the fields by the residual gauge symmetry, the Ward identities impose 
constraints on higher order Green functions due to the same residual gauge symmetry. 

We start defining a Lie operation i^g<i(r,i) by 

V(r,i/(r',t)=:[^'^(r,t),/(r',t)]. (4.5) 

From the algebraic structure of g"-{r,t) ( [4.1| ) we get 

Vi(ri,*)V^(r.,*)^'^^(r3,t) = (if r'^r--^^5\r^ - r^)5\r^ - r^)g\r^,t), (4.6) 

or 

VMri,t)V2(r„t) ■ ■ ■ Vn.i(r„_„i)(7""(r„,t) = F"--'^"'^"+Hri, ■ ■ ■ ,r„)^'^"+nri,t), (4.7) 

where F"i"""""+i(ri, • ■ ■ , rn) is a complex function of {ri, ■ ■ ■ , Fn}. Thus we obtain 
from iWM 

(p|[^»nri,t),^"Hr2,t)]b)=0 (4.8) 

(p|Vi(r„i)V2(r„t)^"^(r3,t)|p) = (4.9) 

and 

(plVi(ri,t)V2(r2,i)--- V"-i(rn-i,t)5'''"(rn,0b) = 0. (4.10) 

The relation between (|4.^ ), ( |4.9| ), ( |4.10|) and gauge invariance can be established as 
follows. Since H due to the temporal gauge has a residual gauge symmetry generated 
by the local algebra {(yf"(r,t)}, the residual group U{g) will generate a set Tigp of 
physical states \p)g from any physical state \p) 

\p)g = U{g)\p)- iov\p)gen,g, (4.11) 

where the residual group element is given by 

U{g) = exp (i f Sr g^uj^\ , (4.12) 

where uo^ = uj^{r) are the local gauge transformation angles. The state \p)g has the 
same energy as \p) since 

g{p\H\p)g = {p\U-\g)HU{g)\p) = {p\H\p). (4.13) 

The weak Gauss law is fulfilled as well 

g{p\g''\p)g = {p\exp{-i J g^u^ £r) g" exp (i J g^ to' £r)\p) 

= {p\Uab9'\p) = 0, (4.14) 

where Uab is a color matrix containing the gauge transformation angles. 
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From the theory of Lie groups it is evident that ( |4.1CI| ) and ( ^.14| ) are equivalent. 



Furthermore, since U{g) generates a gauge transformation for the SU(N) gauge fields, 
any gauge invariant physical observable O has the same expectation value within the 

Since the weak Gauss law is conserved in time within a 4-point truncation scheme, 
all Ward identities which follow from the group properties ( [4.1|) of the Gauss law 
operators in combination with the weak Gauss law, such as ( [4.8|) , (^4.9| ) and (|4.1CI| ), 



will also be conserved. This is due to the fact that correlation dynamics in general 
conserves the commutation and anticommutation relations of all field operators within 
Green functions, which leads to a conservation of the expectation values of all operator 
equations generated by these relations. The problem of the Ward identities thus is 
again shifted to an initial value problem. 

5 Summary and discussion 

In this work we have derived a closed set of equations of motion for connected Green 
functions of quark and SU(N) gluon fields up to the 4-point level. The resulting set 
of equations, which we for abbreviation will denote by Constraint QCD (CQCD) 
equations furtheron, fulfill the conservation of linear momentum, angular-momentum 
and total energy throughout time as well as the weak Gauss law and the Ward identities 
- as generated by commutators of Gauss operators - provided that the initial conditions 
follow the gauge constraints. The latter may be achieved, for example, by a functional 
containing the chromomagnetic field (2.8), which obeys the Gauss law [^. It might be 



used as a trial wave function in a variational calculation within a finite number of basis 
states. Alternatively, one might start from perturbative gluon-field configurations and 
increase the coupling g from zero adiabatically for the preparation of an initial state Q. 

We note that the CQCD equations can easily be transformed into transport equa- 
tions by means of a Wigner transformation with respect to the space-time degrees of 
freedom, which is nothing but a unitary transformation. Since the latter transfor- 
mation has been performed quite often in the literature, we do not present the final 
expressions for reasons of missing compactness. 

The actual goal is a numerical integration of the CQCD equations within a finite ba- 
sis set in analogy to the works performed in the nuclear physics context p^ - |2^ which 



appear to quite naturally explain a variety of non-perturbative nuclear phenomena ^0 
Though the actual expressions in Appendix B appear very cumbersome, we note that 
the reduction of the quark-gluon many-body problem - as e.g. studied by lattice QCD 
in thermal equilibrium - to a nonlinear two-body problem provides a numerical task 
of lower complexity than the original QCD equations on the lattice. Furthermore, the 
CQCD equations allow to study the dynamical evolution of non-equilibrium configu- 
rations of quarks and gluons and to explore the low-energy QCD dynamics, especially 
the linear response of systems close to the ground state. However, it is not clear if the 
actual solutions of our approach will reproduce all the physical phenomena of the full 



^Note that the perturbative vacuum |0} is a physical state in our reahzation of the gauge constraint 
because it fulfiUs (4.4). 
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QCD equations. Before speculating about the convergence of the CQCD equations we 
prefer to analyse the numerical results first. Such work is in progress. 



One of the authors (S. J. W.) likes to thank Prof. U. Mosel for the kind hospitality 
at the University of Giessen where this work was performed. 
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Appendix 

A Cluster expansions 



By taking the functional derivative of the generating functional ( |3.1CI| ) with respect 
to the different sources (cf. Sect. 3) we obtain the relevant cluster expansions of 
Green functions up to sixth order in the gluon, quark, and mixed sector of the SU(N) 
gauge theory. For reasons of compactness we introduce the general particle exchange 
operators 

k k—1 k 

St = (1 + E[^l" + ^2n] + E E ^ln^2,), (A.l) 

n=l n=l j=n+l 

where Vij is the two-body permutation operator. Here i and j represent the set of quan- 
tum numbers characterizing the gauge fields, the corresponding conjugate momenta or 
the quark spinors. 

We note that any application of the Vij operators has to ensure the original order 
of the field operators within the connected Green functions. 

A.l Gluon sector 

In the following the integer numbers between brackets stand for gauge fields A"^ or 
conjugate momenta 11". The cluster expansions then can be written in compact form 
as: 

(1) = (l)c, (A.2) 

(12) = (12),+ (1)(2), (A.3) 



3 

(123) = (123)e + (1 + E ^in)(l)(23), + (1)(2)(3), (A.4) 

n=2 



(1234) = (1234)e + (1 + E ^1") (1) (234), + (1 + E ^2n) (12)e(34), 

n=2 n=3 

+53^1)(2)(34),+ (1)(2)(3)(4), (A.5) 

5 

(12345) = (1 + E^in)(l) (2345), 

n=2 

+5|[(12),(345)e+ (1)(2)(345),+ (12),(3)(4)(5)] 

+ (1 + E ^5n)(l + E P2n)(12),(34),(5) + (1) (2) (3) (4) (5) , (A.6) 

n=l 71=3 
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(123456) = (l + f]^P,„,)[(123)e(456), + (l + Pi4^257^36)(123),(4)(5)(6)] 

n=4 j=l 

+5|[(12)e(3456)e+ (1)(2)(3456),+ (12),(3)(4)(5)(6)] 

6 5 

+^3^(1 + Yl P4n)(l)(2)(34),(56), + (1 + E P6n)5|(12)e(345),(6) 
n=5 n=l 

+ (l + ^P5n)(l + EP2n)(12)e(34),(56),+ (l)(2)(3)(4)(5)(6). (A.7) 

n=l n=3 

The number of terms for the 1-, 2-,..., 5-, and 6-point gluon Green functions consist 
out of 1, 2, 5, 15, 52, and 203 different contributions, respectively. Due to the apphed 
truncation scheme we have neglected the connected 5-point function ( 12345) c in ( |A.6D 
and the terms (1 + En=2^in)(l)(23456)c and (123456)^ in (Q. 



A. 2 Quark sector 

The argument z of a quark spinor u{i) in this subsection stands for the spatial coordi- 
nates and the Dirac, flavor and color indices of a quark i. The cluster expansions then 
can be written as: 

(n(l)^lO) = (..(l)^l'))c, (A.8) 

{u{l)u{2)u{2')u{l')) = {u{l)u{2)u{2')u{l')), + (1 - Vv2')WHl'))M'^H2')),, 

(A.9) 

{u{l)u{2)u{3)u{3')u{2')u{l')) = 

(1 - Vv2' - Vv-y - ^2'3' + V2z{Vvy + Vv2')){u{l)u{l'))M'^>{'^'))M3)u{'i')), 

+ (1 - Vv2' - Vvz'){u{l)u{l'))M'^)u{3)u{'i')u{2')), 

+ (1 - Vv2' - V2'z'){u{2)u{2')),{u{l)u{3)u{3')u{Y)), 

+ (1 - Vvy - V2'^'){u{'i)u{'i'))M^)u{2)u{2')u{l')),. (A.IO) 

The number of terms for the 2-, 4-, and 6-point quark Green functions consist out of 
1, 3, and 16 contributions, respectively. Due to the truncation on the 4-point level we 
have neglected the connected 6-point function {u{l)u{2)u{?))u{?i')u{2')u{l')) c in ( |A.1CI| ). 

A. 3 Mixed sector 

In the mixed sector of the SU(N) gauge theory the Green functions contain quark and 
gluon field operators simultaneously. Here the u{i) represent the quark spinors while 
the numbers 2, 3, 4, and 5 stand for gauge fields or conjugate momenta as in Subsection 
A.l. Since the explicit cluster expansions in the gluon sector are presented in ( |A.2D 



(|A.7|) we have not specified the expansions of the pure gluonic Green functions in 



( [A.ll| ) - ( [A.16| ). The expansions up to the 4-point connected level read: 



{u{l)u{l')2) = {u{l)u{V)2), + (n(lMl'))e(2), (A.ll) 
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(n(l)n(l023) = {uil)uil')23)c + (u(l)n(l'))c(23) + (1 + V23){ui^Hl')2),{3),{A.12) 



{u{l)uil')23A) = (^(1)m(1'))c(234) + (1 + P23 + V24){ui^Ml')2),{3A) 
+ (1 + P24 + P34)(m(1)m(1')23),(4), 



(A.13) 



{u{l)u{l')2345) = {u{l)u{l')),{2345) + (1 + ^ P2„)(^(l)u(l')2)c(345) 

n=3 
5 
+ (1 + E[^2n + P3n] +^24^35) (W(l)w(l')23)c(45), (A.14) 



n=4 



(ii(l)ii(2)ii(2')n(l')3) = {u{l)ui2)ui2')u{l')),{3) 

+ (l-Vv2')Ml)uil')),{ui2)ui2')),{3) + {u{l)u{l')3),{u{'^)ui2')), 

+ {uil)uil'))M^H2')3),], (A.15) 



{u{l)u{2)u{2')u{l')3A) = {u{l)u{2)u{2')u{l'))^{3A) 

+ (1 -n'2')[(^(l)^(l'))c(^(2)^(20)c(34) + («(1)«(1')34),(«(2)m(2'))c 

+ (zx(l)n(l'))c(t^(2)n(2')34), + (1 + P34)(^(l)^(l')3)c(^(2)«(2')4),]. (A.16) 



The original cluster expansions ( [A.ll| ) - ( [A.16| ) consist out of 2, 5, 15, 52, 8, and 
17 terms. Because of the truncation scheme on the 4-point level we have neglected 
{u{l)u{V)23A), in (^^, {u{l)u{l')23A5)c and (1 + En=3^2n)(t^(l)t^(l')345)e(2) in 
{\KJ^ , {u{l)u{2)u{2')u{l')3)c in i\KJ^ , and {u{l)u{2)u{2')u{l')3A) ^ and 
{1 +Vs^){u{l)ui2)u{2')u{l')3),{4:) in {Wm. 



B Equations of motion for connected Green func- 
tions 

B.l Gluon sector 

Whereas the time evolution of the 1-point and 2-point gluon Green functions is given 
in ( p.20D - ( p.24| ) we here continue with the 3- and 4-point functions. In order to allow 
for a compact presentation we again use the permutation operators Vn' and Ted as 
defined in Section 3. The equations of motion for the connected 3-point gluon Green 
functions read: 



j^{Atix)At:ix')At::ix")). = (n«(x)A^;(x')A^;;(x"))c 
+(A^(x)n^;(x')A^:(x"))c + {A^ix)A^ix')ir^:ix")), 



(B.l^ 
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f^{mx)At, ix')A1„ (x"))c = (nr(x)n^, ix')At„ (x"))c + {m^)A'l, ix')Ut„ ix'% 



+d'AAt{x)A1, {x')A1„ (x"))c - d,^d^Mti^)At' {x')A1„ (x"))c 

+2gr'^dy,[{Alix)A^iy)At:ix')Ai:ix")), 

+(1 + Z,){{Aiix)At:ix')At:ix"MAKy))c 

+(1 + %,a''){Al{x)At:{x'MAt{y)A:}''W'))c}h=y) 

+g\r"'r'\i + %b + Te,)[(At(x)Af (x)A^;(x')A^;;(x")).(A^(x)), 

+(1 + Tu){Al{x)At:{x')),{Ai{x)Ai:{x"))Ml{x)). 

+{Ai{x))Mi{x))M^{x)Ai{x')Ai:{x")), 

+(1 + T,,,„){A\{x)Ai{x)Ai:{x'))Mt{x)Ai:{x")), 
HA\{x)At{x))M1{x)At{x')Ai:{x"))^] 

+9r"'dy,[{Al{y)A^{x)Ai{x')Ai::{x")), 

+ (1 + %,){{Al{y)AUx')Ai:{x"))M-{x))c 

+ (1 + T„,„.)(4(2/)^-''(^'))c(A-(a:)A«;;(x"))c}](.=,) 

-gr"'dy\{Al{x)Al{y)AUx')4'{x"))c 

+(1 + T,,)(4(a;)yi^;(x')A^;;(x"))c(A^,(y))c 

+(1 + %,a"){A\{x)Ai{x'))Ml{y)Ai:{x")),]\^=y) 



+gt''"'i\pK{x)ul{x)At, {x')At. ix")) 



(B.2) 



£ 
It 



a' / ^l\-na" / tl\ 



(jit{x)iit, {x')At„ (x"))c = {i\-nmt' {x')nt,. {x")) 



+(1 + Vu'WMnx)^' {x')At. (x"))c - d^^d^Ml^mt' (^')^-" (^")) 

+2^(1 + ru')r"'dy,u\{x)j^Xy)^i{x')^''{x"))c 

+(1 + %,){{A\{x)iii,{x')A!t::{x"))M'i{y))c 

+(1 + T„,„.)«(^)nr(x'))c(A^(i/)Ar;(x"))c}](.=,) 



x[{A\{x)Ai{x)YYi:{x')Ai::{x"))Ml{x))c 

+{i+rM){A\{x)uUx'))Mi{x)Ai:{x"))Mnx))c 
+{A\{x))Mi{x))c{Ai{x)ni{x')At::{x")), 

+(1 + T,,a"){Al{x)Ai{x)Ui,{x'))Mt{x)Ai:{x")). 
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+{Ai{x)Ai{x))MK^)m'i^')A^:'i^"))c] 

+9(1 + Vu')r'%A{Aliy)At{x)Ut:{x')At:ix")), 

+ (1 + Z,){{Al{y)Ut:{x')At:{x"))MK^))c 

+(1 + %'a''){Ai{y)U^:{x'MAKx)A'i::{x''))c}]i.=y) 

-g{l+Vu')r'%A{Al{x)Al{y)Ut:{x')Af,{x")), 

+(1 + Zb){{Ai{x)ui:{x')Ai:{x"))Miiy))c 

+(1 + T,,,.){Al{x)W}:{x'))Ml{y)Af,:{x")),}]^,=y) 



+^(1 + V,,^t^"'i\^{u'^{x)u%{x)JYl, {x')A% (x"))e, (B.3) 

^(n«(x)n^;(x')n^:(x"))c = 

(1 + v., + p,,0[5^(^- (^)n^''(^')n?:(x"))c - a.^a.,(A^(x)n^;(x')n^;;(x"))c] 



+(1 + T,,){(yit(x)n«;(x')n^;;(x"))c(A^(?/))c 

+ (1 + T,,,,0(4(^)n?'(a:'))c(A^(?/)n^;;(x"))c}](.=,) 
x[(4(x)Af(x)n^;(x')n^;;(x"))c(A^(x))e 

+(4(^))c(A^(x)).(A^(x)n^;(x')n^;;(x"))c 

+(1 + T,,')i,A\{x)Ai{x)Xli{x'))M\{x)Ai:,{x")), 
-rm^)Ai{x))^{A\{x)W(,{x')W(:,{x"))^\ 

+(?(! + P,, + V.,.)r"'dy,U\{y)A\{x)W:,[x')J\%{x")), 

+(1 + T,,){(4(y)n^;(x')n»;;(x"))e(A^(:r)), 

+(1 + T„,,„)(4(i/)nr''(^'))c(A-(x)n«;;(x"))c}](.=,) 

+(1 + T,,){(4(x)n^;(x')n«;;(x"))c(A^(i/))c 

+(1 + T.,..)(4(^)nr''(:^'))c(A^fc(?/)nr:(x"))c}](.=.) 

+(7(1 + V,v + P..")^"''7a/3(«a(a:)«M^)nr''(a^')nr'' (a^"))c. (B.4) 

The equations of motion for connected 4-point gluon Green functions are the most 
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cumbersome of our approach and read: 

+(An^)^"'(^')^?'(^")n«;;;(x'"))c, (b.s) 



c 



+2gr'^dy,{i + T,,)[(4(x)A^;(x')A^;;(x")A^;;;(x'"))c(A^(?/))c 

+ (1 + T,,„,. + T„,v'0(4(^)^"'(^0^"''(^''))c(^^(?/)Ar;"(a:'''))c](.=,) 

x[((<(x))e(A^(x))e + (<(x)Af(x))e)(A^(x)A^;(xOA^;;(x'OA^;;;(a;''o 
+(1 + T,,){(At(x)A^;(xO).(A^(x)A^;;(a:"))e(A^(:r)A^;;;(x'"))c 

+ (1 + T,,,,. + T„,w'0(4(^)^"''(2:0^"''(2:'0)c(^'(x)A^:'(x'''))c(A'(x))e} 



c 



+(1 + T,,,,„ + T,,v'0(4(^)A''(^0A''(^''))c(A(?/)A''' (^'''))J(.=.), (B.6) 



I (n^(x)n^; (x') a: (x") a^;;; (x'"))c = (n^(x)n^; (xon^: (x'o a;;; (x-)). 

+ (n»(a;)n»;(a;')A;;(a;")n^;;;(x'"))c 



+ (1 + V.,)[dt{At{x)Iit> (x')A„ (a;")A"' (^'"))c 
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A{{Al{x))Mi{^))c+{Al{x)Ai{x))^){Al{x)uUx')Ai:^ 

+(1 + T,,){(4(a;)n^;(x')).(Af (x)A^;;(x"))c(A^(x)A«;;;(x'"))c 
+(1 + T,,,. + T„,„.o{(^^(^)^f (^)nr''(a:'))c(A^(x)A^;;(x")A^;;;(x'"))c 

+(4(^)^^(^)^?'(^'0^-"'(^'''))c(An^)nr''(^'))c}] 



+^(1 + p..)r'^5..(i + %,)ui{y)Y\^:{x')Ai:{x")Ai'::{x"'))M'i{^))c 
+(1 + T,,,,. + T„.„.o(4(2/)n^'\2:o^-"(2:''))c(A'(x)A^;;;(x'''))d(.=,) 

-i7(i + ^..or^ii + rMAi{x)iiUx')Ai:{x'')Ai::{x'''))Mi{y))c 

+ (1 + T,,,,. + T„.„.0(4(2^)n^''(2^0A"(2:''))c(A(?/)A'" (x'''))J(.=,), (B.7) 



-(n^(x)n^, (x')n:'„ {x")At,„ (x'"))e = (nr(x)n^, (x')n^. (x")n^., (x'"))c 

-a.,9.,(A^(x)n^;(x')n^;;(x")A^;;;(x'"))c] 

+2^(1 + vw + n.or^^9,,(i + T,,)[(4(a;)n«;(a;')n»;;(x")A^;;;(x'"))c(A^(?/)). 



x[((4(^))c(A^(x))e+(4(^)^f(:^))c)(A-(x)n^;(xon^;;(a:'OA^;;;(x'''))e 

+(1 + T,,){«(x)n^;(x'))c(Af (a:)n^;;(x"))c(A^(x)A^;;;(x'"))c 

+(1 + T„v" + T„.,„o(4(^)n?'(x')n^;;(x"))e(Af (x)A^;;;(x'")).(A^(x)) J 

+(4(x)A^(a;)n^;;(x")Ar;;;(a:'"))c(A^(x)n^;(x'))c}] 

+<7(i + v., + p.,or''=5..(i + r,,)[(4(y)n^;(x')n^;;(x")A^;;;(x"'))c(A^(:r)), 

-(?(! + n, + ru")r"'dyXi + r,,)[(4(x)n^;(a;')n^;;(x")Ar;;;(x'"))c(A^(i/))c 

+(1 + %,,,. + T„.,,„)(4(^)n?'(a;')n^;;(a;"))c(A^(l/)A^:'(:r'"))J(.=,), (B. 



|(n«(x)n^;(x')n^;;(x")n^:;(x'"))c 



:i + Vw + Pn" + Vu'")mA1{x)Yit, {x')nt. {x")ni„ {x'")) 
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a' ( u\ rra" ( ^11 \ rra" 



-d.^d^Mii^mt' {x')iit„ {x")Yit,. ix'")) 



+2gil + Vw + Vu" + Pu"')r'%,a + Zb)[{Al{x)Uf,{x')Uf,{x'')U^:{x'''))Mtiy))c 

+(1 + %,a"' + T„.,.0(4(^)n?''(^')n?'' (^"))c(^^(i/)nr;;;(x'"))c](.=,) 

x[((4(^))c(Af(:r)),+ (4(x)Af(x))J(A^(a;)n^;(x0n^;;(a:'0n^:'(x'''))c 

+(1 + r,,){(4(x)n^;(x'))c(A^(a;)n^;;(x"))c(A^(x)n^;;;(x'")). 

+(1 + %.,.. + T,„,,„)(^^(^)nr''(^')n?" (^"))c(Af (a;)n«;;;(x'"))e(A^(:r)) J 
+(1 + T,.a" + T,,,.0{(4(^)4(^)nr''(^'))c(4^(x)n^:(a:")n^:(:r'")). 

+g{i + p., + n.. + Vu'.)r''dy,{i + r,,)[(4(2/)n^;(x')n^;;(x")n^:'(a;'"))c(^-(^))e 
+(1 + %,,,. + T,.,.0(4(?/)nr'(x')n^:(x"))c(A'(^)n^r(x'"))c](.=,) 

+(1 + Ta.a"' + T,.,.o«(a^)nr'(^')n?'''(^"))c(^^.(y)nr:'(x"'))c](.=,). (b.9) 

B.2 Mixed sector 

The remaining equations of motion - not specified in Section 3 - read as follows: 

i{ul{x)ui,{x')Yi\{y)4{y')), + i{ul{x)ui{x')A\{y)YiUy'))c 

+<K(^^(^)<'(^')^-(?/)4'(?/'))c 
-9t'''^^m{x)ui,{x'))M':{^)A\{y)AUy'))c+{ul{x^^^^^^^ 

+al,p[d^>^{ul{x)4{x')A\{y)4{y')), 

-gr'''^{{ul{x)ul{x'))Mi:{^')A\{y)AUy'))c + {<{x)^^^^^^ 

+ {<{x)u%x')4{y'))M':{x')A\{y))^ 

+ «(x)«^(x')4'(y)4'(?/'))c(<(x'))e}], (B.IO) 

^J^{ul{^)<.{x')4{y)AUy'))c = i{ui{x)ui,{x')4{y)4{y'))c 
+4J,d^,{u%{x)ui{x')n\{y)AUy'))c 

~gr^^{{u%x)ui{x'))M':{^n{y)A%')).+ {ul{^^^^^^ 
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-^t^'''^'^{(<(x)«^(xO)e(<(xOn^(y)4\?/0)e + (<(x)«^(xOn^(y)).(<(x04'(y')). 

+ (<(x)«^(xOA^:(yO)c(<(xOn,^(y))c+(<(a:)^z^(xOn,^(y)A^:(yO)c(<(a:'))c}] 



+i2gr^d-.,[{ul{x)ui{x')A':{y))M-{^)Almc 

+ (<(x)<:(xOA^(.)).(<(l/)4'(l/'))e+(<(x)<:(x')A^(^)4'(?/'))e(<(2/))c 

+«(x)<:(x')4"(i/)4'(?/'))c(4^(^))c]fe=.) 

+tg'f'"^r\l + T,., + T,,g[«(x)<:(x')4"(2/)).(4(y)A^(?/)4'(?/'))c 

+ (l+^.e)«(x)<:(x')4"(?/))c(Af(y))e(^-(2/)4'(?/'))c 
+ (l+T,.,T,,e)«(x)<:(x')4"(|/)4'(l/'))e(^f(2/)^?(?/))c 
+ «(x)<:(x')4"(?/)4'(2/'))e(Af(y))e(^^(y))e] 

+Z^f^'X[«(^)<''(^')4"(^))e(4^(2/)4'(2/'))e 

+ (<(x)<:(xOA^(?y))e(<(^)4'(l/'))e+(<(x)<:(x')4^(2/)4'(2/'))e(4'>))c 

+ «(x)<:(x')4'»4'(?/'))c(4^(?/))c](,=.) 

-^^r'X[«(^)<'(^')4"(2/))c(^^.(^)4'(?/'))e 

+ «(x)<:(x')A^(^))e(4"(2/)4'(l/'))e + «(x)<:(x')A^(z)4'(l/'))e(4"(l/)). 

+ «(x)<:(y)4"(?/)4'(2/'))e(A^(^))c](,=.) 

+z^t^"'^7«';3"^(^' - y){ul{x)ul,{y)AUy'))c 

+ «(^)«^'(l/))e(4"(l/)<:(x')4'(?/0)c], 



(B.ll) 



-gr'^^{{ui{x)ui,{x'))M':{m\{y)^Uy'))c + {u%^^^^^ 
+(«^(x)<:(xOn^:(2/0)c(4'\^)n^(y))e + (t^Kx)<:(xOn,^(2/)n^:(i/'))e(4''(:^))c}] 

+«''/.[5.',«(^)<(^')n,^(2/)n^'(?/0)c 
~gr'^'^{{ui{x)ui{x'))Mi:{^'n{y)^Uy'))c + {^^^^^^ 
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+(<(x)«^(xOn^:(yO)c(<(xOn,^(y))c + (<(a:)«^(x')n,^(i/)n^:(?/'))c(<(:r'))e}] 

+z(l + Vu')\dl{ul{x)ui,{x')A\{y)Ii%'))^ - dy^dy,{ul{x)ui,{x')A\{y)IiUy'))c] 

+(<(x)<:(x')A^(^))e(4''(?/)nr(2/'))e + (<(x)<:(x')A^(^)n^:(y'))c(<(?/))c 

+ «(x)<:(a:')4"(?/)nr(2/'))e(4^(^))c](,=.) 

+^9\l + ^..)r"7^'^'^(i + %., + r,,g[«(x)<:(x')4"(?/))c(4(2/)^^(?/)nr(2/'))e 
+(i + r,e)«(x)<:(x')4"(?/))c(Af(i/))e(A^(2/)nr(?/0)c 

+ (l+T,.,T,,e)«(x)<:(x')4"(2/)n^''(l/'))e(^f(l/)^^(2/))c 
+ «(x)<:(x')4"(?/)n^:(2/'))c(^f(2/))e(^^(l/))e] 

+^^(l + P..')r'X[«(^)<'(^')4"(^))c(^■(?/)^^:(2/'))e 

+(<(x)<:(xOA^(?/))e(4'>)nr(i/'))e+(<(x)<:(x')4^(?/)n^:(?/'))e(4'>))c 

+«(x)<:(x')4"(^)n^'(?/'))c(4^(2/))c](,=.) 

-^^(i + ^..')f''X[«(^)<''(^')4"(?/))c(A^(^)nr(i/'))e 

+«(x)<:(y)4"(2/)nr(i/'))e(A^.(^))4=.) 

+ «(^)«^'(l/))e(4"(l/)<:(x')nr(?/'))c]. (B.12) 

This completes the nonhnear CQCD equations, which form a closed set of first order 
differential equations in time and can be integrated by standard numerical techniques 
for almost arbitrary initial conditions, provided that these fulfill the weak Gauss law 
as well as the Ward identities as discussed in Section 4. 
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